An unusual space-time evolution for heavy ion collisions at high energies 

due to the QCD phase transition 
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The space-time evolution of high energy non-central heavy ion collisions is studied with relativistic 
hydrodynamics. The results are very sensitive to the Equation of State (EoS). For an EoS with the 
QCD phase transition, an unusual matter distribution develops. Before freeze-out, two shells are 
formed which physically separate and leave a maximum in the center. We make specific predictions 
for the azimuthal dependence of the flow and for two-pion interferometry, contrasting our results 
with a resonance gas EoS. 



1. One of the principal goals of the heavy ion colli- 
sion program is to find and to quantify the QCD phase 
transition from hadronic matter to a new phase, the 
quark-gluon plasma (QGP) Q. Experiments at the 
Brookhaven AGS (lab energy 11 A*GeV) and at the 
CERN SPS (lab energy 200 A*GeV) are expected to 
produce a QGP/mixed phase during the initial stages 
of the collision, although currently there is only indi- 
rect evidence for this state (see the recent reviews |H). 
With the completion of the Relativistic Heavy Ion Col- 
lider(RHIC) at Brookhaven and its much higher collision 
energy (100+100 GeV*A in the center of mass frame), 
the experiments are expected to produce the QGP well 
above the transition temperature. In this work, we study 
how the strong QGP pressure can be observed at RHIC. 

The position-momentum correlations of the produced 
hadrons, colloquially known as collective flow, directly 
reflect the EoS of the excited matter. Multiple studies 
usingcascade event generators and hydrodynamics (see 
e.g. [|]-H )) have successfully reproduced the AGS/SPS 
hadronic spectra. A radial flow velocity of about (0.5- 
0.6)c is found in central PbPb collisions but the flow 
develops principally during the late hadronic stages of 
the collisions and has little to do with the QGP. An EoS 
extracted from these model studies shows "softness" dur- 
ing the early stages of the collision, either due to the 
proximity of the QCD phase transition J5j, or due to 
non-equilibrium phenomena such as the formation and 
fragmentation of strings || . 

2. Additional information about the EoS may be ex- 
tracted from the azimuthal dependence of flow in non- 
central collisions, which depends non-trivially on the im- 
pact parameter and the collision energy. The ellipticity of 
the flow has been studied theoretically [|7|-|lCfl and exper- 
imentally pHJl^| . Because elliptic flow develops earlier 
than radial flow, its systematic measurement at the SPS 
may settle the mixed phase/pre-equilibrium controversy 
mentioned above. The original purpose of the present 
study was to further quantify elliptic flow within a hydro- 
dynamic framework. Instead, we found that non-central 
collisions at RHIC/LHC energies have an unusual expan- 
sion pattern, which cannot be described as simply elliptic 
and which is qualitatively different from AGS/SPS ener- 
gies. 



3. Let us begin with a description of the transverse 
acceleration. In model calculations, the radial accel- 
eration history changes from SPS to RHIC due to the 
QCD phase transition p3[ . The ratio of pressure to en- 
ergy density, p/e, has a deep minimum at the end of 
the mixed phase, known as the "softest point" of the 
equation of state JI3J]. For AGS/SPS collision energies 
the matter is produced close to the softest point and the 
resulting transverse acceleration is small. Therefore, in 
non-central collisions the matter retains its initial elliptic 
shape and burns slowly inward. For RHIC /LHC collision 
energies, the early pressure starts an outward expansion 
|fl3f . This outward expansion and the inward deflagration 
can cancel each other, making a stationary front, called 
the "burning log" in p5[ . Summarizing, at the AGS/SPS 
there is first softness and then a hadronic push, while at 
RHIC/LHC there is first a quark-gluon push, then soft- 
ness, and then a hadronic push. In spite of this change, 
the final radial flow velocities at RHIC and at the SPS 
are expected to be similar JTgJ . 

4. The early push redistributes the matter, however. 
The early velocity has long time, ~ 10/to/c, to influence 
the matter distribution before freeze-out. The stiff QGP 
in the center, with T >> T c , pushes against the soft mat- 
ter on the exterior, with T m T c , producing a shell-like 
structure. Since the final distorted distribution rather 
resembles a nut and its shell, we call this picturesque 
configuration the nutshell. For high energy non-central 
collisions, the matter expands preferentially in the im- 
pact parameter direction (the x axis) and the expanding 
shells leave a rarefaction behind. Furthermore, since the 
acceleration started rather early, the two half-shells par- 
tially separate, and by freeze-out three distinct fireballs 
are actually produced. We have called this consequence 
of early pressure the nutcracker scenario. 

5. Following we assume a rapidity-independent 
longitudinal expansion. We then solve the 2+1 dimen- 
sional relativistic Euler equations in the transverse plane, 
with the coordinates x,y and proper time r = (t 2 — z 2 ) 1 / 2 , 
using the HLLE Gudunov method |18| |. As in previous 
calculations flTif , we have used a simple bag model equa- 
tion of state with T c = 160M eV and a lGeV latent heat. 
We have modeled hadronic matter with a simple reso- 
nance gas EoS, p — .2e [O]. The pressure was taken 
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FIG. 1. Typical matter distributions in the transverse 
plane at mid rapidity, calculated using boost invariant hydro- 
dynamics, for PbPb collisions at b=8 fm and proper time r in 
a representative RHIC collision. For the corresponding cen- 
tral collision the pion multiplicity, including all isospin states, 
was dN^/dy = 850. The outer dotted line and inner dashed 
line show temperatures of 120 and 140 MeV respectively. The 
solid lines show contours of constant energy density with a 
step size of 10 MeV. (a) shows the distribution at RHIC for 
a resonance gas EoS, p = .2e. (b) shows the distribution at 
RHIC for a bag model EoS. 

to be independent of baryon number which is a good 
approximation at high energies. The initial entropy dis- 
tribution in the transverse plane was assumed to be pro- 
portional to distribution of participating nucleons as in 
0. We parameterize the initial energy density by the 
total pion multiplicity, dN^/dy. (How particle multiplic- 
ity maps to the collisions energy depends on the entropy 
production mechanism. This mapping will soon be deter- 
mined experimentally at RHIC.) For dcfmitcness, we con- 
sider PbPb collisions at an impact parameter of b = 8/m 
and freeze-out at a fixed temperature, Tj — UOMeV. 

At SPS energies the flow develops late, and the matter 
retains its initial almond shape until the late hadronic 
stage. However, at RHIC the flow develops early and 
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FIG. 2. The distribution (in arbitrary units) of the x coor- 
dinates of nucleons emitted from the spatial band, \y\ < 2/m, 
and integrated over intervals of proper time, (a) shows the dis- 
tribution for a resonance gas EoS. Grey, hatched and dark his- 
tograms are for r < 8, 8 < r < 9 and 9 < r < (T max = 10.5) 
fm/c , respectively, (b) shows the distribution for a bag model 
EoS. Grey, hatched and dark histograms are for r < 10, 
10 < r < 11 and 11 < r < (T ma x ~ 12.5) fm/c, respectively. 

redistributes the matter by late the hadronic stage. Two 
sample matter distributions, in the transverse plane and 
at fixed proper time, are shown in Fig. [|. A resonance 
gas EoS, p — .2e, produces little structure and simple 
elliptic flow (Fig. |l]a). An ideal gas EoS, p = e/3, also 
produces simple elliptic flow and even shorter lifetimes 
(not shown). Finally, our bag model EoS produces two 
"nut-shells" of matter which expand outward (Fig. [j]b). 
Note that in Fig. [l]b the matter is pushed into two shells, 
moving in x direction, while at the north and south poles, 
two holes develop. A maximum, the "nut", remains in 
the center. The matter distribution becomes increasingly 
"nutty" with larger impact parameters and higher colli- 
sion energies. 

The evolution is clarified by plotting the emission 
points of nucleons, integrated over periods of proper time. 
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In Fig.g, (a) and (b), we show the x coordinates of the 
emitted nucleons for a resonance gas EoS and for a bag 
model EoS respectively. For both EoS, during the early 
stages, particles are slowly emitted from from a station- 
ary freeze-out surface, making a peak at x = 4 — 5 fm 
||15| . For both EoS, however, 75% of the nucleons freeze- 
out during a short proper time interval of 2.5 fm/c. For a 
resonance gas EoS (Fig. ||a) the freeze-out positions are 
uniform and become increasingly centralized with time. 
In contrast, for a bag model EoS (Fig. ||b) the distri- 
bution has three distinct and comparable sources for the 
final 1.5 fm/c. There is a central "nut" and two extremal 
"shells" with a hole in between. 

6. We turn now to the experimental consequences of 
this flow pattern. First we examine the <f> distribution of 
the produced particles at mid-rapidity (y=0). These dis- 
tributions are expanded in harmonics and are sometimes 
weighted by the transverse momentum squared. 
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We have calculated the single particle distributions for 
various secondaries, using the standard Cooper-Frye for- 
mula f2(i|j . The elliptic components, v-2 and depend 
only weakly on collision energy, as found in previous stud- 
ies [Q. For nucleons, for example, we found v 2 « 7% and 
a 2 ~ 13% from the highest SPS energies to LHC energies. 

Higher harmonics, in contrast, grow from SPS to 
RHIC. To summarize the effects of higher harmonics in 
the distributions, we have plotted in Fig. || the weighted 
net nucleon <j> distribution (the l.h.s. of equation(2)) for 
SPS and RHIC, normalized to the first two terms in the 
Fourier expansions shown above. In the dashed curve cor- 
responding to RHIC, the early pressure forces a 3 — 4% 
additional asymmetry in final net nucleon distributions 
beyond the elliptic component. At LHC energies the addi- 
tional asymmetry is even more pronounced. The marked 
minimum at 45° is due to the square shape of the mater 
distribution, which somewhat reduces the flow on the 
diagonal. Note also a prominent positive correction to 
elliptic flow at 90°. Both of these effects are observable, 
given the expected statistics at RHIC. 

The distribution of deuterons and other heavy frag- 
ments should express the underlying flow more clearly. 
Because the emission points of the nucleons are bunched 
along the ridges of the nutshell, larger fragments are gen- 
erally emitted from the shells. This inhomogeneity en- 
hances the production probability of fragments and peaks 
their final flow in the x direction. Multiply strange barons 
such as f2~ are also of interest. Since they do not re- 
scatter in the hadron phase, they reflect the early flow 
pl| . Indeed any azimuthal dependence of the flow of mul- 
tiply strange baryons would be fairly convincing evidence 
of collective motion in the quark phase. 
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FIG. 3. The net nucleon <j> distribution weighted by 
(the l.h.s. of equation(2)) and divided by the first 
two moments of the corresponding Fourier expansion (the 
r.h.s. of equation(2)). This plot summarizes the effect of 
all higher moments. The solid curve is for the SPS and 
the dashed curve is for RHIC. In the corresponding cen- 
tral SPS and RHIC collisions the pion multiplicity, including 
all isospin states, was dN^/dy — 525 and dN^/dy = 850 
respectively. We summarize the numerical values below. 
At the SPS : < pi >- 67 GeV 2 , a 2 = 12.0%, a 4 = 0.6% 



At RHIC 
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<Pt >-- 



GeV 2 
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Q2 
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Note that ao is proportional to < p\ >. 



7. The spatial asymmetry of the matter distribution 
at freeze-out is probed by Hanbury Brown- Twiss (HBT) 
two particle interferometry. Strong flows strongly mod- 
ify the source function. Each correlator with given pair 
momenta, is generated by its own "patch" , or "homo- 
geneity region" |22| . Taking these patches together gives 
a complete picture of the source. We will discuss this 
complicated issue elsewhere, and here show only two se- 
lected correlators which emphasize the qualitatively dif- 
ferent predictions of different EoS. 

The correlators are found by taking the appropriate 
Fourier transform of the source function over the freeze- 
out surface ^0|. Below we display the Hanbury Brown- 
Twiss(HBT) radii R xx and R yy and employ the notation 
(pi +P2)ti — 2A" p and (pi -p 2 ) l i = ftu- For R xx we select 
gin the x direction , the direction we want to probe. K is 
chosen in the orthogonal y direction , with magnitude .5 
GeV. For R yy the axes of K and q are simply reversed. 
The correlators may then be fit to the functional form 
C = 1 + exp^qfRft), where Rf^ has been interpreted as 
the source size at zero velocity 123]. More specifically 
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These radii are shown in Fig. ^ for PbPb collisions at 
b=8 fm, as a function of the pion multiplicity scaled by 
the number of participants to central collisions (not b=8 
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FIG. 4. HBT radii, R xx and R yy , for \K\ = .5GeK, as a 
function of the total pion multiplicity. The quoted multiplic- 
ity includes all isospin states and is scaled by the number of 
participants from b=8 fm to b=0 fm (a factor of 2.9) . See 
the text for more details. For a resonance gas EoS p = .2e, 
the solid lines and dotted lines show R xx and R yy respec- 
tively. For a bag model EoS, the dashed and dashed-dotted 
lines show R xx and R yy respectively. 

fm). R xx and R yy are shown for a bag model EoS and 
for a simple resonance gas EoS, p — .2e. For low energies 
near the left hand side of the plot, the two EoS show 
approximately the same radii, roughly corresponding to 
the initial elliptic shape of the matter distribution. For 
a simple resonance gas EoS, the HBT Radii show lit- 
tle energy dependence, while for an EoS with the phase 
transition the homogeneity regions increase steadily with 
beam energy. The rapid increase of R xx can be under- 
stood qualitatively. The contributing pions move in the 
y direction with rather high momenta, .5 GeV. The pair 
therefore originates, not from ridges of the nutshell, but 
from the region in between the shells. The rapid increase 
of R xx reflects the increased separation of the nutshells 
at higher collision energies. The increase in R yy reflects 
the flattening of the shells themselves. For flat, square- 
shaped, shells the homogeneity regions are larger than 
for curved elliptic shells. 

8. In conclusion, for non-central heavy ion collisions, 
we predict an unusual space time evolution, which re- 
sults from the interplay of a hard and soft EoS typical 
of the QCD phase transition. We have called this the 
"nutcracker" flow since two shells are produced and then 
separate. The azimuthal momentum asymmetry can be 
seen in the higher harmonics and flow of heavy secon- 
daries, while the spatial asymmetry can be seen in HBT 
interferometry. 

We end with the experimental strategy. As the 
"nutcracker" flow persists for all sufficiently non-central 
events, and because the principal RHIC detectors can de- 



termine the impact parameter plane, absolutely any ob- 
servable, from strangeness, to flow, to J/ip suppression, 
should display marked azimuthal dependence, which re- 
flects the fireball in different stages. We urge our exper- 
imental colleagues to look for this dependence from the 
first day of operation at RHIC. 
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